Exact Bond Ordered Ground State for the Transition Between the Band and the Mott Insulator 
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We derive an effective Hamiltonian Hef / for an ionic Hubbard chain, valid for t <C f/, A, where t is the 
hopping, U the Coulomb repulsion, and A the charge transfer energy. Heff is the minimal model for describing 
the transition from the band insulator (BI) (A — f/ S> t) and the Mott insulator (MI) (f/ — A » t). Using 
spin-particle transformations (Phys. Rev. Lett. 86, 1082 (2001)), we map HeffiU = A) into an SU(3) 
antiferromagnetic Heisenberg model whose exact ground state is known. In this way, we show rigorously that a 
spontaneously dimerized insulating ferroelectric phase appears in the transition region between the BI and MI. 

PACS numbers: 71.27.+a, 71.28.+d, 77.80.-e 



It is well known that the insulating state can have differ- 
ent origins. The simplest case is the BI since it is one of the 
possible solutions for non-interacting electrons moving in a 
periodic potential. In particular, the BI occurs when the num- 
ber of particles per unit cell is even, so the bands are either 
full or empty. The other traditional example of insulating 
state is the MI. In this case, the charge localization is just a 
consequence of the local Coulomb interaction U between the 
electrons. The charge gap of the BI is just the band gap. In 
the case of the MI, this gap is a function of the Coulomb re- 
pulsion and goes asymptotically to U in the strong coupling 
limit. These two insulating states have completely different 
properties. For instance, the BI's are paramagnets while the 
Mi's, in general, are antiferromagnets. This observation rises 
the following question: what happens when a system evolves 
continuously from the BI to the MI state? 

To find an answer to this question has been the main mo- 
tivation for studying the Ionic Hubbard model (IHM) during 
the last ten years. The IHM was originally proposed in the 
80's to describe the neutral-ionic transition in mixed stack 
charge-transfer organic crystals 01 Ql' This model is a Hub- 
bard Hamiltonian on a bipartite lattice with different diagonal 
energies for the two sublattices. The difference between both 
energies is A. During the 90's, there was a renewed interest 
in the IHM due to the potential applications to the description 
of the ferroelectric (FE) perovskites QHIH. At half filling, 
the ground state of the IHM is an ionic or BI for A ^ J7 and 
a MI for [/ ^ A. Field theory arguments in one dimension 
pointed out the existence of an intermediate bond ordered in- 
sulating (BOI) phase for A < C/ < i 0]- For i < C/, A, 
perturbation theory clearly describes the BI (i ^ A — J7) and 
the MI {t <^ [/ — A) including the charge dynamics of 
the latter However, perturbation theory diverges in the 
transition region and no insight is provided for the BOI. The 
numerical solutions of finite chains have difficul- 

ties and contradictory conclusions were reported. The reason 
will become clear after deriving the main result of the present 
paper. 

Contrary to the cases of the BI and the MI states, no ex- 
act solution having long range bond ordering is known for the 
transition regime. Finding an exact ground state is not only 
crucial to prove the existence of the BOI phase rigorously. 



but also to understand its microscopic origin and fundamen- 
tal properties. This result becomes even more important if we 
consider that the BOI phase is an electronically induced spin- 
Peierls instability that generates a FE state out of the spin- 
singlet dimer pairs. In particular, a bond ordered FE state was 
observed in the pressure-temperature phase diagram of the 
prototype compound, tetrathiafulvalene-p-chloranil fil l fl2il . 
In addition, as it was pointed out by Egami et al. the 
microscopic origin of the FE transition in covalent perovskite 
oxides like BaTiOa is still unclear. It is known that a pic- 
ture based on static Coulomb interactions and the sirnple shell 
model is inadequate to describe some FE properties I13ll . The 
exact result presented here demonstrates that when an ionic 
insulator gets close to a charge transfer instability (strong co- 
valency), an electronic mechanism for ferroelectricity takes 
place. 

In this paper, we first derive an effective Hamiltonian, 
Heff, for the limit [/ > f and A > i of an extended IHM that 
includes a nearest neighbor Coulomb repulsion V . By means 
of the generahzed spin-particle transformations introduced in 
Refs Ullillil, we map Hef / into an anisotropic SU(3) an- 
tiferromagnetic Heisenberg model that becomes isotropic for 
U = IS. and particular values of the ot her p arameters. The 
isotropic model is exactly solvable flTl Il8l fl^ and is also 
equivalent to the biquadratic S — \ Heisenberg model. The 
exact ground state is a dimerized spin system that becomes a 
BOI when translated back to the original fermionic variables. 

We will start considering an IHM with an additional 
Coulomb interaction V between nearest neighbors: 

i,(T i 

+ U J2{n^^ - i)(na - ^) 
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+ yE(n,-l)(n,+i-l), (I) 
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where = fj^fi^, = J^a ^i'^^ t > 0. Note that A > 
is the diagonal energy difference between sites in different 
sublattices. 

We will only consider the half-filled case p — I, i.e., one 
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particle per site. If t = 0, the ground state of if is a non- 
degenerate BI for U < A + 2V. There is only one low-energy 
state per site: the odd sites are doubly occupied and the even 
sites are empty (see Fig[0. For U > A + 2V the ground state 
is a degenerate MI (one particle per site) and the low-energy 
subspace has two states per site due to the spin degeneracy. 
When t is finite and small, far from the transition region the 
low-energy Hamiltonian of the MI is a Heisenberg model and 
the degeneracy is lifted in favor of a spin-density wave 01 
(see FigOJ. However, to describe the transition region, we 
need to include three states per site in the low-energy subspace 
Hq. In order to construct an effective Hamiltonian in Ho, it is 
convenient to perform an electron-hole transformation for the 
odd sites; 
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After this transformation, H becomes invariant under a trans- 
lation of one lattice space {i ^ i + 1). The low energy 
subspace TYq is now defined as the set of states with no dou- 
ble occupancy on any site. This constraint can be incorpo- 
rated by defining the constrained fermion operators: = 
cL(l — cLc.^f), and c- = fcL)^. In addition, to connect our 
low energy theory with spin Hamiltonians (see below), we in- 
troduce the following transformation: 



^4 



for i = An and z = 4?! + 1 



t for i = 472 + 2 and i = 4n + 1 (3) 



The high-energy subspace 7i i consists of all states which have 
at least one double occupied site. The only term of R that 
mixes Tip and TYi is the hopping term t. By means of a sec- 
ond order canonical transformation that eliminates the terms 
proportional to t, we obtain the following effective Hamilto- 
nian for the low energy spectrum of H: 



i 
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y^(l-fi,)(l-n^+i), 



(4) 



where rii = J^a 4a(^'ia and sf = 1/2 Y,r,T' At'^tt'<^^t' with 
a = {x, y, z] (cr" are the Pauli matrices). The exchange in- 
teraction J = 2t^/ ([/ + A — t^) comes from a second order 
process in the hopping t. In this derivation we have neglected 
the second order three-site hopping term. The negative sign 
for the XX and yy interactions is just a consequence of the 
gauge transformation of Eq. (|3}. 

\it ^ V and A — U ^ t, the ground state of H^f / is the 
empty state that corresponds to the BI in the original language. 
For [/ — A ^ t, the ground state has one particle per site. 
By eliminating states with empty sites with another canonical 
transformation, H^ff reduces to the Heisenberg model that 



FIG. 1: Schematic plot of the different ground states of H, 



eff- 



To exploit the symmetries of Heff, it is convenient to use 
the generalized Jordan- Wigner transformations introduced in 
Refs. Ellillil. In particular, for this case it is appropriate 
rewrite Hef / in terms of the hierarchical SU(3) language l!l6|l. 
To this end, we need to map the constrained fermion operators 
in one sublattice (say A) into SU(3) spins in the fundamental 
or "quark" representation JTslfl^ : 



4^ 



"jT - 3 



where is the kink operator H 

Ki = exp[i7r 



k<j 



nk 



(5) 



(6) 



that transmutes the statistics. The components S^'^ are 
generators of the su(3) algebra with commutation relations 

= V'^'^^^'O') - V'5''^'(J)- We also 
make use of the conjugated or "anti-quark" representation: 



-J "it ^iiST '"-jT - 

to describe the degrees of freedom of the B sublattice. Eqs. (|5} 
and (0 are generalizations of the Jordan- Wigner transforma- 
tions to SU(3) spins The resulting H^ff is as an 
anisotropic SU(3) Heisenberg model with an applied "mag- 
netic field": 



(7) 
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with Joo = -V, Joi = Jo2 = -t, Jii = Ji2 = J22 = 
-J/2, J^^ J^^, and B = 2^/3 - J/3 + U/2 - A/2. Note 
that this model connects SU(3) spins in the A sublattice with 
the conjugate SU(3) spins (anti-quark representation) in the B 
sublattice. Let us now consider V ~ J/2 = t and U — A. For 
this line in the space of parameters Hef / is an isotropic SU(3) 
antiferromagnetic Heisenberg model that is invariant under 
staggered conjugate SU(3) rotations, TZ and TZ\ on sublattices 
A and B respectively. This model is integrable ITtI ITM fl9ll 
and the exact solution is a spin-dimerized ground state. In the 
original language this is equivalent to saying that the charge 
and the spin are both dimerized (bond ordering). The exact 
ground state energy per site is e^/t = —1.796864.... The 
value of the gap is A/i = 0.173178, a rather small value, 
and the correlation length ^ = 21.0728505... is very large 
fl^ . This explains the numerical difficulties for identifying 
this phase. 

To make clear the relation between the dimerization in 
SU(3) and the bond ordering of the fermionic variables, we 
just need to translate the corresponding order parameter from 
the SU(3) language back to our original fermionic language. 
The spin-dimer SU(3) order parameter is I.21i,i22i1 : 

i5= (9) 

where = |i — 1, i) (i, i — 1| is a projector on the SU(3) 

singlet spin state, \i — l,i), at the bond (t — 1, j). Note, in 
addition, that /li-i.i is the isotropic SU(3) Heisenberg Hamil- 
tonian for the same bond. In the fermionic language, |i — 1, i) 
has the following expression: 

In terms of the original /-fermions [see Eqs.(|3 and (|3}], this 
is a linear combination of an on-site singlet and a nearest- 
neighbor singlet state that is illustrated at the bottom of Fig. 
Xian showed that D = 0.4216IJn where Dq is the value of 
D for a perfect dimerized state 1I21I1 . Replacing ( I10> in (|9}, it 
becomes clear that D is a BOI order parameter From a calcu- 
lation of the charge Berry phase |7] we obtain that the shift 
in polarization of the two possible perfect dimerized states 
(see the bottom of Fig^ relative to the MI is ±e/6. Eq. (I10> 
shows that the dimer formation is just a consequence of the 
charge transfer instability between the two sublattices. In this 
sense, this exact ground state unveils the fundamental role 
of covalency for the stabilization of a bond ordered FE state 
This is not only relevant to describe one dimen- 
sional systems like tetrathiafulvalene-p-chloranil jTllfT2il . but 
also helpful to gain a deeper understanding of the FE transi- 
tion of covalent perovskites ||3|. 

Using another set of transformations that connect the con- 
strained fermions with SU(2) 5=1 spins lll4ll : 

5+ = x/2(ct^X,+if]c^.,), 

SJ = V2iK}c^,+cl^K,), 

Sj = ^jT -riji, (11) 

we can also write Hf,ff{V = J/2 = i) as an S = 1 bi- 
quadratic Heisenberg model fisll with a single-ion anisotropy: 



H^ffiV = J/2 = t) = -t^ (S, . S,+i)' + i?^(5f )^ 

i i 

(12) 

where E ^ {A—U)/2. The strength of the anisotropy term is 
determined by the difference A — U. If E' is large and positive 
(A ^ U), the spin system has an easy plane anisotropy. Each 
site is most of the time in the S*^ = state, which means that 
the magnetization is perpendicular to the z-axis. The ground 
state is non-degenerate (there is no broken symmetry) and cor- 
responds, in the original language, to the band insulator. If E 
is large and negative [U ^ A), the system has a strong easy 
axis anisotropy and each site is in the — ±1 state, i.e., 
the local magnetization is parallel to the z-axis. The ground 
state is critical due to the antiferromagnetic correlations which 
characterize the MI. In between, for U = A, we have demon- 
strated that there is a dimerized state which corresponds to the 
BOI. In terms of the original variables, the strong quantum 
fluctuations that appear in the proximity of the charge transfer 
instability break the Z2 inversion symmetry by increasing the 
strength of one bond relative to the next one. 

The S — 1 version of H^ff also provides a simple way 
of studying the low energy excitations of our exact bond or- 
dered state. The excitations of a dimerized spin 1/2 chain are 
spinous that carry a spin S = 1/2. In the same way, the 
excitations of our S = 1 dimer state are 5=1 spinous. 
Each spinon is a soliton or anti-phase boundary for the Z2 
spin-dimer order parameter The two regions with opposite 
dimerization are separated by a local 5=1 defect which is 
attached to the anti-phase boundary. In terms of the origi- 
nal language, the 5z = spinous correspond to s = soli- 
tons (charge excitations), while the 5^ = ±1 spinous are 
s = 1/2 solitons (spin and charge excitations). From Eqs. 
(|3 and il2\ . the total charge operator relative to half-filling is 
Q = 'EtGAi^!)^ - Y^i^Bi^tf- For a general dimerized so- 
lution with arbitrary E we have on each site: ((5f )^) = 1 — a 
with Q = 1/3 for the exact solution al E = {). When the 
defect is localized on site j (the extension of the defect does 
not affect the charge or spin of the excitations because they 
are topological invariants) we have ((5|)^) = 1,0. Since 
the defect spearates two regions with opposite dimerization, it 
is easy to check that the s = solitons (anti-solitons) have 
charge Q = ±(1 — a) while the charge of the s = 1/2 
solitons(anti-solitons) is Q = ±a. For E = 0, both excita- 
tions are degenerate as a consequence of the SU(3) invariance. 
These excitations coincide with those obtained by Fabrizio et 
al 01, who treat the bosonized IHM as a phenomenological 
Ginzburg-Landau energy functional. The magnitud a is pro- 
portional to their jump in the charge field. In particular, the 
s = 1/2 excitations interpolate between an ordinary electron 
near the BI-BOI boundary (a = 1, £■ > 0) and a spinon near 
the the BOI-MI boundary {a = 0, E < 0). 

In order to obtain an exact ground state of H^f / we used 
a value of J = 2t which is beyond the region allowed by 
perturbation theory. Therefore, to connect our exact solution 
with the IHM, it is important to discuss the effect of reduc- 
ing the value of J. In addition, since most of the previous 
papers do not include the V term, it is also important to an- 
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alyze its effect. A simple first order estimation of the energy 
change for the three competing phases, SEbi = SU /2 + t, 
SEmi = 2tln2, SEboi = - 5eot/9, when J and V 

are reduced to zero, indicates that the BOI survives if the dif- 
ference U — Ais simultaneously increased to a value of order 
t. This conclusion is supported by different numerical results 
QUI^ that report a bond ordered ground state of H{V = 0) 
for the same region of parameters {U — A ^ t and U,A^ t). 
These observations suggest that our exact solution is contin- 
uously connected with the BOI phase which was numerically 
found in the IHM. 

In summary, we have derived an effective low energy 
Hamiltonian, i?e//, for the t/, A > i limit of the IHM. Hef / 
is a minimal model to describe the BI to MI transition. Its 
simple form and the fact that it operates in a reduced Hilbert 
space of local dimension D — 3 provide a new framework to 
understand this transition. Using the spin-particle transforma- 
tions introduced in Refs. fisl fl6ll . we mapped Hefj into an 
anisotropic SU(3) antiferromagnetic Heisenberg model. By 
increasing the value of J beyond the region allowed by pertur- 
bation theory, we have shown that there is an exactly solvable 
SU(3) invariant point for [/ = A and V = J/2 = t. In this 
way, we demonstrated the existence of a bond ordered phase 
for the transition regime between the BI and the MI. The large 
value of the correlation length ^ = 21.0728505... explains 
the numerical difficulties for detecting this phase in finite size 
systems. This exact solution provides an answer to the ques- 
tion that was formulated in the introduction: in the transition 



regime the system resolves the competition between the BI 
and the MI by creating a rather local resonance which can be 
visualized on each dimer as a linear combination of a unit cell 
of the BI and a nearest-neighbor singlet that is related with 
the MI (see Fig^. One of the most important physical conse- 
quences of this non-trivial phenomenon is the generation of a 
new mechanism for ferroelectricity f^k that should be relevant 
for ionic insulators which are close to a charge-transfer insta- 
bility. This result should motivate a careful reexamination of 
the ferroelectricity in covalent materials. 

In addition to the SU(3) mapping, we used a second 
spin-particle transformation lll4ll that maps the constrained 
fermions into 5 = 1 SU(2) spins. In this new language 
Heff{V = J/2 = is a biquadratic Heisenberg model with 
a single-ion anisotropy along the z-axis which is proportional 
to A — L/. The transitions between the BI, BOI, and MI phases 
have been re-interpreted in the S ~ 1 language. We have also 
shown that the S = 1 version of He/ / provides the most natu- 
ral frame to understand the low energy excitations of the BOI. 
These excitations are solitons that carry a fractional charge 
and spin s = or s = 1/2. The relevance of these results 
illustrates the elegance and the potential of the generalized 
spin-particle transformations introduced in Refs. lll4lll5Lri6ll . 
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